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Abstract 

We calculate the phase space volume occupied by a nonextensive system of ./V classical 
particles described by an equilibrium (or steady-state, or long-term stationary state of a nonequi- 
librium system) distribution function, which slightly deviates from Maxwell-Boltzmann (MB) 
distribution in the high energy tail. We explicitly require that the number of accessible microstates 
does not change respect to the extensive MB case. We also derive, within a classical scheme, an 
analytical expression of the elementary cell that can be seen as a macrocell, different from the 
third power of Planck constant. Thermodynamic quantities like entropy, chemical potential and 
free energy of a classical ideal gas, depending on elementary cell, are evaluated. Considering the 
fractional deviation from MB distribution we can deduce a physical meaning of the nonextensive 
parameter q of the Tsallis nonextensive thermostatistics in terms of particle correlation func- 
tions (valid at least in the case, discussed in this work, of small deviations from MB standard case). 
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I. INTRODUCTION 



Statistical description of a system of iV particles requires the subdivision of the phase 
space into equidimensional elementary cells of phase volume A Q, which can be determined 
by the laws of nature (comparison with quantum evaluation of energy state density) and 
experimentally measured, for instance, in the low temperature heat capacity of a crystal or 
in the Stefan-Boltzmann constant. 

In the phase space volume of a system of particles described by quantum distribution the 
smallest elementary cell is the third power of the Planck constant. For classical particles 
the elementary cell is, in principle, undetermined. This is true particularly in the limiting 
case of small occupation numbers (when MB distribution is valid) and the phase space 
volume of a cell acquires arbitrary values. 

The problem of dividing the phase space into finite cells was solved when a natural way 
to derive dimension of elementary cells within MB statistics using energy quantization 
was found. However, quantum discontinuity can be lost when dimension of the volume 
containing the gas increases and the quantum states become even and even more numerous 
[lj]. Only the insertion of the Pauli exclusion principle enables us to solve this problem. 



Boltzmann's request of a great number of partic 
particle every 30000 cells in normal conditions) 



es into the cell is usually not verified (one 
However, if the number is small, it is 



possible to combine many cells and form a greater cell (macrocell) containing more particles 



3, 



This paper aims at examing first of all how , in the nonextensive thermostatistics (NETS), 
the elementary cell differs from the one of the extensive MB case, requiring explicitly that 
the number of accessible microstates be the same in both (extensive and nonextensive) phase 
spaces and obtain also, within a classical scheme and without quantum arguments, explicit 
expressions of the cells. Let us briefly. recall that NETS has been developed in the last years 



eny n 



mainly after the works of Tsallis j^, |5[ which introduced a generalised entropy, featured by 
an entropic parameter q, whose meaning is not fully understood yet (more generally, it can 
be defined in terms of fluctuations of intensive quantities like temperature jy, Q, ; explicit 
analytical expressions of a in terms of physical quantities exist in few cases, like, for instance, 



in turbulence problems [9j and in the description of solar plasma 



10 



111). 
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The value of q is not easily given a priori for a given particular system, without comparing 
calculations to available experimental results. 

The NETS is now applied to a great variety of problems, from genetics and biology to 



astrophysics (see ref.s 



121 and 



13j for a complete list of basic works and applications). 
The approach to NETS we are showing in this paper, is based on the analysis of deviations 
from standard phase space volume and on a new definition of elementary cell. This new 
approach will probabily provide, in the near future, a better understanding of the meaning 
of the parameter q. The validity of our treatment can be extended to statistics different from 
the Tsallis version of NETS, when complete deformed distribution functions are considered, 
and to large deviations from MB phase space volume. 

In Section 2, by comparing the number of allowed microstates Wo in the MB extensive phase 
space Qq and in the deformed phase space Q, we derive the dimensions of both elementary 
cells AQ and AQ. The cell AQ results to be smaller than AQ (if q < 1) so that, in this 
case, the third power of Planck constant is not suitable as the value of Af^o because h 3 is 
the smallest cell admissible due to the Heisenberg principle. Otherwise AQ is larger than, 
AQq if Q > 1 {h is expressed in units [energy length]). 

By considering the nonextensive Tsallis equilibrium (or steady state, or long-term stationary 
state of a nonequilibrium system) distribution, we limit ourselves to small deviations, with 
the absolute value of the parameter 5 = (1 — g)/2 much smaller than one (we recall that for 
q — > 1 all the MB results must be recovered). 

We must realize that deviations, although small, are not negligible in those energy islands of 
Q where quantities like, for instance, nuclear reaction rates, rates of atomic processes and of 



chemical reactions, electron transport in semiconductors are very sensible to them jiol 

In Section 3, we discuss the case of nonextensive classical ideal gas. We com par e 



14. 



18. 



our results with those recently derived by several authors, among them see ref.s |l7 . 
and we calculate the entropy, the chemical potential and free energy discussing the results. 
We deduce also in Section 4 an interpretation of the parameter q using the calculated phase 
volume Q as a function of fractional deviation, occupied by a nonextensive system of N 
classical particles. The physical meaning is derived in terms of particle correlation function. 
Conclusions are reported in Section 5. 
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II. PHASE SPACE VOLUME AND ELEMENTARY CELL 



The state of a system of N particles is specified in the 6 N dimensional phase space (r 
space). Its volume Q[n r ] can be divided into small cells of volume AQ, so that coordinates 
do not vary sensibly within them. 

The volume contains the N particles distributed in a certain set of numbers 
(ni, n 2 , ■ ■ • , n t ) — [n r ] and is given by: 

n[n r ] = —— (AQ) N . (2.1) 

Hi! n 2 \ ■■■ , n t \ 

The number of accessible microstates is defined (including the factor 1/iV!) by 



n\ (Any 



W = . (2.2) 



It is well known that In Q can be represented through the Stirling approximation 

InQ^i N \nN + N InAQ-Y.nilnrii . (2.3) 

i 

Having Q a maximum, its logarithm can be expanded in a series of powers of drii = rij — n io , 
where n io is the distribution function for which Q has a maximum, indicated by f2 (the 
symbol d means variation). After using the Lagrange method with the usual constraints 
and with negligible interactions the following is obtained: 

3 

n ki = A M e~ Xi , n *o = N > Hi o x * = 9 N > 

i i 1 

x i = (3e l , /3= YT' (2 ' 4) 

V \27rm J 

The mass of the particles m is in energy units. The number of accessible microstates in Qq 
is of course 

w ° = «T(f?vr • ' 2 ' 5 ' 

(let consider the following example: for a classical ideal gas, if iV is large we have that 



W = e iV /v / 2WV, see Sect. 3). 

We are now interested in the volume in the T space corresponding to a NETS distribution 
having small deviations from the MB distribution. The deviations in the high energy tail 



are particularly interesting. Following the approach to general statistical problem in physics 



by Bohm and Schiitzer 



21 1 after expanding up to the second power in drii the Lagrange 



equation, alternative definition of Q is obtained: 

lnfi = lnfi -i 

° 2 V ru 



(2.6) 



or 



Q = Qq exp 



n io ( d fi. 



(2.7) 



with the fractional deviation from MB distribution given by 

drii 



(2.8) 



Let us introduce the non-Maxwellian distribution 



n-, 



A T A s exp (-yi - 5 yfj 



(2.9) 



Eq. (j2.9j) represents a distribution, mainly differing from Maxwellian in the high energy tail, 
depending on the sign of 5 [5 > 0, (q < 1): depleted tail; 5 < 0, (q > 1): enhanced tail]; it 
can be derived, for instance, from the Tsallis distribution 



n T = Z q (P 



l-(l-q)P" (E^-^ 



1-9 



(2.10) 



when deviations from MB distribution are small, Z q {j3*) is the generalized partition function 
and U q the internal energy. Eq. ()2.9|) may represent an equilibrium distribution or a steady 
state, or a long-term stationary state of a nonequilibrium system. In Eq. (12.91) 

3/2 



.4, 



N ( P* 



An 



(2.11) 



V \2Tim / 

AQ is the deformed elementary cell of the phase space Q, f3* is the Lagrange multiplier of 
the phase space Q, different from (3 (the Lagrange multiplier in Qq), defined below in ()2.24|) . 
We also have 

15 



A 5 = 1 + ^5 - 30 5 2 

Vi = P* €i ■ 



(2.12) 
(2.13) 
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£j are the energies of the different microstates with the same values of those of the phase 
space Qq. The constraints are 

3 



E 



ri; 



N 



N 



(2.14) 



where the second one is the energy NETS average value when deviations are small. Other 



22 



2311 : of course, when deformations are 



distributions, based on other statistics can be used 
small, these distributions can be fitted with good approximation by Eq. (j2.9|) . The use of 
exact distributions, derived from other generalized statistics is possible at least numerically; 
in this work, however, we limit ourselves to consider only the distribution of Eq. ()2.9j) . in 
order to simplify the question. 

The volumes Qq and Q are different in size. AQ and AQ are also different. We realize 
that the differences could be considered negligible, while their effects seem to be quite 
important for the evaluation of several physical quantities. Let us remark that the number 
of microstates or discrete events does not change from one space to the other one. Therefore, 
we set the equation W = Wo, because we want to count the same number of microstates 
both in Q and Qq. Using the relations reported above in this Section , we can explicitate 
the expressions of Q and AQ after simple calculations. After using the following relations: 

AQ ( 8*y e 95 / 4 



Ai 



AQ 



A 2 s 
£3 



1 



n? 



o 

1185 
N I I 



5 2 



n 



Hi 



18 5 



A s 



3799 
~32~ 



5< 



(2.15) 
(2.16) 
(2.17) 



noting that the explicit calculation of ()2.17|) requires the evaluation of the average value 

/e" Xl \ /37 c 11841 r2 \ 

= « Hi'— sr*) • (2 - 18) 

finally, we obtain Q and AQ as functions of Qq and Aficb respectively as follows: 

1 



Q 

Q~ ~ : 

AQ 

~AQq~ 



exp 



-In 

2 



N 



E 



exp 



1 



3799 



exp 



2 V 32 



,18 5- 
2 V 32 



5< 



(2.19) 
(2.20) 



The systems featured by A > (q < 1) have AQ smaller than A^o- 

Therefore, we are not allowed in their classical extensive description to take AQq = h 3 
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because h 3 should be the smallest elementary permissible cell. 

On the other hand, we shall verify that AQq must always be much larger than h 3 . It must 
be a macrocell. 

The equation of state calculated by means of the distribution (|2.9|) is given by 

PV = NkTCg = NC s /(3 , (2.21) 

where 

C 5 = 1 -55 + AQ5 2 . (2.22) 

Let us anticipate that for a nonextensive classical ideal gas after exact calculations with the 



exact distribution (|2.1(Jj) . one can obtain the following relation, at any value of q 171 Il8] 



PV = Nf3* = NC g /P , (2.23) 
P*C q = P, (2.24) 

where C q = Z q ((3y- q . 

We compare ((HHP and fFT^i to (El and (fiZty in the case of small deviations. With 
the condition 3iV(l — q)/2 = 3N5 > 1 the quantity C q reduces to C$ if we assume an 
elementary cell given by 

and we recover the correct formal expression of the equation of state reported in Eq. (j2.21|) . 
This value of AQ imposes that 

3 1 221 

A u = l, and A T = l--5+—5 2 . (2.26) 

Therefore, the requirements that W = W and that the equation of state for NETS classical 
systems be correctly expressed also in the small deviations limit imply that the standard 
phase space elementary cell be given by the expression (j2.25j) (we send to Sect. 3 for some 
more details). 

This elementary cell can indeed be considered a macrocell, particularly if compared to the 
value of h 3 , usually taken as elementary cell. However, this requirement is not a problem, 
because of the uncertainty of the classical elementary cell and because Darwin and Fowler jjj 



7 



^2.27) 



showed that macrocells should be used to satisfy Boltzmann requirements of a great average 
number of particles in each cell. 

To take AQ Q of Eq. (|2.25|) means to have one particle in each macrocell and, posing AQ = 
Xh 3 , to have 1/X particles in each microcell h 3 , i.e. one particle in thousands of cells, 
where 

/27rm\ 3/2 V_ 1_ 
Nh?' 

which is not a pure number, but depends on (3. 
The elementary cell in the deformed phase space is from ()2.2()j) 

which is a macrocell. We can write AQ = Yh 3 , each macrocell containing one of the N 
particles, where 



^2.28) 



Y 



'2vrm\ 3/2 V Cf 2 



(3 ) N h 3 ' 
is a quantity depending on (3 and q or on (3*. 

The subdivision of the phase space volume Q in microcells equal to h 3 does not allow the 
conservation of the value of the number of permissible microstates, nor the subdivision in a 
fixed number of microcells independent on (3 and q. 

Instead of having an elementary cell which does not depend on the Lagrange multiplier f3 and 
is the third power of a universal constant (Planck constant) and the number of elementary 
cells forming the total phase space volume depending on the Lagrange multiplier (3, we have 
an elementary cell which depends on the Lagrange multiplier (3 (it is a macrocell). We also 
have the number of cells depending only on the number of particles iV (with these positions 
we may satisfy the requirements to conserve the number of accessible microstates both in 
extensive and in nonextensive phase volumes Qq and Q and to preserve the correct form of 
the equation of state). 

III. APPLICATION TO CLASSICAL IDEAL GAS: ENTROPY, CHEMICAL PO- 
TENTIAL, FREE ENERGY. 



Classical ideal gas model based on nonextensive thermostatistic,^ 
subject of several studies since the first applications of NETS \24 . 



relations has been the 
25^ . Classical ideal gas is 
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described by an unperturbed state of a system with long-range interaction and the model 

n 

can be solved analytically [171 ] . In the past the equation of state has been derived within 
several schemes along the steps of evolution of NETS: non-normalized |24j, |25j], normalized 
17 1 , OLM (optimal Lagrange multiplier) incomplete statistics 0,15 among oth- 
ers. 

In this Section we show synthetically that the formally correct equation of state ()2.23p ob- 
tained by means of the exact expression of the distribution ()2.10|) . can be deduced also in 
the small deviations case, by taking into account the distribution ()2.9|) if the macrocells of 
Eq. ()2.25j) and of Eq. ()2.28)1 are taken as elementary cells. 

The equation of state specific for classical ideal gas remains form invariant under nonex- 
tensive generalization of thermodynamics. That is, it is valid for all q: PV = N//3*. This 
means that the equation of state of a NETS classical ideal gas is formally equal to that of 
extensive statistical mechanics: P V = N/ (3. 

The Lagrange multiplier f3* previously introduced, associated with the constraint in NETS, 
is defined by 

/?* = ir4— > 

-t phys 

&t is a constant depending on q which becomes the Boltzmann constant k for q — > 1 |3| 
and Tphys, the physical temperature, is 

Tphys = w = = i 1 + hf s <) k^p ' (3 - 2) 

(Sj is the Tsallis entropy defined below in Eq. ()3.1()j) : in the first treatment of classical ideal 
gas based on OLM approach [3] the dependence on q was attributed only to and not to 
^phys or P*)- The partition function Z q (j3) is defined by 

V N r 

Z M = TTTTXTTw / II d *P* n T , (3.3) 



where is the square bracket factor of Tsallis distribution [Eq. (j2.1(Jj) ]. 

Two other quantities that are useful for NETS classical ideal gas calculations are: the 

coefficient 

V N r n q 

c « = ]«WwWn^^ = W-«, (3.4) 
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and the internal energy 



(3.5) 



1 V f 3 

^ = C£ A! (Afio)^ 7 ^ d3p< ? 



2 m Z q {fiY 



From the explicit expression of Z g , C 9 and C/ g (we do not report it here and we send the 
reader to ref.s we ma y verify that the above three functions Z q , C q and U q depend 

(when q ^ 1) on the elementary cell AQq as 

zm « (Aa,)-^ 1 -^ , 

C^iAQo)^^ 1 -^ 1 -^ , (3.6) 
C7g ps (Afi ) _Ar(1 ~ 9)/(1_Q) , 

where Q = 3JV(1 -q)/2. 

As well known, when q —>■ 1 we have that Zi ~ (Afio) _7V an d C\ and C/i do not depend 
on AQq- The above functions enter into the calculation of the equation of state, which can 
be derived by means of the usual thermodynamic relations. As we have already discussed 
in the previous Section, by considering a nonextensive classical ideal gas with distribution 
function rii of Eq. ()2.9j) (small deviations from MB distribution), we have calculated that 
the equation of state is given by PV = N kT (1 — 5 5 + 46 5 2 ) i.e., in the limit of small 
deviations we must have (3* = /3/C$ as it can be easily verified. 

In fact, within the treatment illustrated in the previous Section, we obtain that in the limit 
of small deviations and for 3 N 5 > 1 the expression of C q reduces to C$ if the elementary 
cell (macrocell) Af2 of Eq. (|2.25j) is assumed. 

With the expression of AQ given by (|2.28J) in place of AQq given by (|2.25J1 we obtain that 
Z q , C q and U q do not explicitly depend on the elementary cell. Z q and C q do not depend on 
P either: 



rffef) n n ( i \ 3N/2 



l + (l-g)|jV 



1-9 1 2 



N 



(3.7) 



and C q and U q can be calculated from the relations 



(the functions T(x) can be calculated by means of the relation T(x) = \J2 tt x x ~ 1 ^ 2 e~ x ). 
Let us note from Eq. (j2.1|) that space phase volumes fl and fl depend on the elementary 

10 



cells Afl and Afl, respectively. If AQ is a constant, like h 3 , Q does not depend on /3, but 
only on N and Q depends on N and q. Instead, if AQq has the expression of Eq. ()2.25|) and 
AQ is given by Eq. (|2.28|) . then fl is also a function of (3 and Q is also a function of f3 and 
Q- 

Let us now derive one of the thermodynamics quantities depending on AQ, the entropy. 
The Boltzmann entropy is defined by 



S B = k In Wo 



- N k + const. 



- N k + N k In 
2 



V 


(2nm s 


3/2 1 


N 


V , 





(3.9) 



where (2 n m/ f3) 3 / 2 / AQ , is the single particle number of cells. 

Taking the elementary cell AQq given by Eq. ()2.25|) . we obtain Sb = 5Nk/2 because the 
single particle number of cells equals N and the value of the constant is zero. 
The nonextensive entropy is 



r-T 



where the condition W 



q 



(w 1 - 



l-q 

k 



(wt q - 1 

Sb 



exp 



k 



l-q) 



- 1 



Wn has been used and therefore 



5 Nk + (l-q)^N 2 k = SB (l + lsN 



(3.10) 



(3.11) 



q 2 

This relation shows that Sj tends to change as iV 2 with large N. The chemical potential \x 
of a classical ideal gas is an intensive quantity, defined by 

N AQ 



/i = kT In 



(3.12) 



V (2nmkTyi 2 _ 

where AVIq is the elementary cell taken usually equal to h 3 . The quantity \x can be negative 
or positive and is equal to zero only at the particular value of temperature 

-iV\ 2 / 3 h 2 



kT 



V 



2 Tim 



(3.13) 



Considering an ideal classical gas, if we increase the number N of one unity adding, at 
constant energy and volume, one particle with kinetic energy zero and neglecting all interac- 
tions, the number of accessible microstates increases and entropy increases too. Considering 
the relation 



AU = T AS — P AV + \i AN 
11 



(3.14) 



we obtain that \i must be negative. On the contrary, if repulsive interactions are active, the 
incoming particle increases both the energy of the system and the entropy; the system must 
spend some energy to return to the initial energy value. Then 5b decreases and sometimes 
decreases more than the gained quantity because of one more particle in the system. In this 
case \i is positive. 

Let us now take for AQq the definition of Eq. ()2.25|) . We find that for an extensive ideal 
classical gas we have fi = 0. This result can be explained because adding one particle at 
constant energy and volume the elementary cells decrease their single volumes and the total 
work done by the system to diminish the spatial volume of N + 1 cells is k T. Therefore, 
from Eq. (j3.14|) we have 

fj,= -kT --kT + kT = . (3.15) 

When the ideal classical gas is nonextensive, using (J3.ll)) the expression of [i is 

45 25 

fis = -—5kT-Y$NkT + 466 2 kT, (3.16) 

which gives \x = when q — > 1 (5 — 0) and is slightly positive or negative depending on 
the sign of 5. In this case the quantity fi is composed by intensive and extensive terms. Its 
behaviour versus k T depends on 5 and on 5 N which can be finite or can approach zero or 
infinity depending on 5 and N, separately. To complete our study let us discuss the free 
energy F. For an extensive ideal classical gas with elementary cell AQ given by Eq. ()2.25|) 
at constant energy and total volume, the function F and its variation due to the addition 
of one particle to the system are given by the following relations: 

F = U — T S B = —N k T , (3.17) 
AF = —kT . (3.18) 

being \i = and considering the work spent by the system to change the spatial volume of 
elementary cells. 

The variation of free energy of a nonextensive ideal classical gas due to the addition of one 
particle to the system is, using ()3.1(i|) 

/ 25 25 \ 

AF q = -PAV + fi s = -kTC s + fi8 = -kT (1 + —5+—5N) . (3.19) 
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This expression differs from that of the non normalized treatment of Ref. 
free energy (as also defined in Abe et al. [18] and with Sj of Eq. (j3.10J) is 



. 3- 



The NETS 



Fq = Uq ~k^ = ~ NkT ( 1 + T 5iV ) ' (3 - 20) 

and AF q equals the same quantity calculated above [Eq. (J3.19)) ]. A new definition of free 
energy introduced by Abe et al. |l8( which is 

F q = U q -T^lnC q = U q -T^ln(l + ^-S q ^ , (3.21) 

gives, by using the elementary cell we have proposed, the same result obtained in Eq. (|3.2(J|) 
and then AF q of Eq. (j3~P3jh 



IV. INTERPRETATION OF THE PARAMETER q 

Let us consider again the square fractional deviation from a MB distribution 

= 2™ f) N , (4.1) 

i rij V 32 / 

We have already discussed that appreciable fractional deviations are, in general, very small. 
Nevertheless, few islands of phase space can give a non negligible contribution to deviation of 
quantities calculated within the standard space, like, for instance, rates of nuclear reactions 
or atomic processes or rates of chemical reactions. We define, as usual, the radial correlation 
function 

9(r) = x J2 6 ( ¥ + T i- ¥ j) > ( 4 - 2 ) 

(g(r) = 1: perfect gas). 
By using the relation 

E^ = l + ^ f\s(T)-l]dr, (4.3) 

from Eq. (|4.1|) we can write the following relation: 

iQ l-q 3799 (1 -qf 1 1 f. . x , lA AS 

18 — = 1 / \q(r) - 1 ] dr . (4.4) 

2 32 4 N V J K 1 

Disregarding, for simplicity, the term in (1 — q) 2 we have 

q = 1 ~l {jf + V J \9(r)-l]dr] . (4.5) 
13 



Therefore, q — > 1 if both the conditions g(r) — > 1 and A?" 
N — > oo but the system is not a perfect gas [g(r) 7^ 1]: 



00 are fulfilled. Otherwise, if 



q = 1 ~lv I [ ,/(r) 



11 dr 



(4.6) 



and if the system is a perfect gas (g = 1) but the number of particles is small and finite: 

1 



Q 



9N 



(4.7) 



Finally, in the limit of a gas/fluid of N particles, we have found a physical interpretation of 
the nonextensive parameter q in terms of the particle correlation function among others (see 



Wilk et al 



3, Beck Q, 



Beck and Cohen 



and Tsallis 29]). It can be considered valid at 



least for values of q not too different from unity, (small deviations from MB distribution). 



V. CONCLUSIONS 



We have calculated the phase space volume f2 and the corresponding elementary cell 
AQ, occupied by N classical, identical particles composing a weakly nonextensive system, 
at equilibrium (or in a steady state, or long-term stationary state of a nonequilibrium 
system). Both quantities must differ from their extensive MB values Qq and AQq because 
we require that the number Wq of available microstates in the MB extensive phase space 
Qq equals the number W in the deformed (nonextensive) space: Wq = W. This condition 
imposes particular analytical expressions of AQq and AQ. The number of elementary cells 
results equal to the number of particles N. 

The two elementary cells result to be macrocells if compared to the usually taken elementary 

cell volume, equal to the third power of the Planck constant , with consequences on the 

expression of thermodynamic quantities depending on elementary cell, like, for instance, 

the partition function, constant of entropy, chemical potential and free energy. 

We have applied results on elementary cell to nonextensive classical gas, recovering known 

results (equation of state) in NETS literature where AQ is usually taken equal to h 3 or left 

undetermined. 

We have derived expressions of entropy, chemical potential and free energy and discussed 
more deeply the case of chemical potential. Applications to particular systems and physical 
cases like nuclear fusion reactions in plasmas and chemical reactions will be reported 
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elsewhere. Through expressions of Q and AQ we have calculated fractional deviation from 
Maxwell-Boltzmann distribution. A physical meaning of the nonextensive Tsallis parameter 
q in terms of the correlation function has been derived. 

We wish to tank Dr. A.M. Scarfone for critical reading of the manuscript and comments. 
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